Wigner crystal in snaked nanochannels: outlook 
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Abstract 

We study properties of Wigner crystal in snaked nanochannels and show that they are characterized by a conducting 
sliding phase at low charge densities and an insulating pinned phase emerging above a certain critical charge 
density. We trace parallels between this model problem and the Little suggestion for electron transport in organic 
molecules. We also show that in presence of periodic potential inside the snaked channel the sliding phase exists 
only inside a certain window of electron densities that has similarities with a pressure dependence of conductivity 
in organic conductors. Our studies show emergence of dynamical glassy phase in a purely periodic potential in 
absence of any disorder that can explain enormously slow variations of resistivity in organic conductors. Finally we 
discuss the KAM concept of superfluidity induced by repulsive Coulomb interaction between electrons. We argue 
that the transition from the sliding KAM phase to the pinned Aubry phase corresponds to the superfluid-insulator 
transition. 
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1. Introduction 

The Wigner crystal [1] appears when the energy of 
Coulomb repulsion between charges of same sign be- 
comes dominant comparing to kinetic energy of charge 
motion. On a one-dimensional (ID) straight line this 
crystal can move ballistically as a whole at an arbi- 
trary small velocity. Here, we discuss the properties of 
Wigner crystal sliding in Id in presence of a periodic 
potential and in Id snaked nanochannel following re- 
cent works [2,3]. An example of snaked nanochannel 
of a sinusoidal form is shown in Fig. 1. The snaked 
form of a channel is very similar to the Little sugges- 
tion [4] on possibilities of superconductivity in organic 
molecules. As described in [4,5], it is assumed that or- 
ganic molecules form some effective wiggled or snaked 
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channel with an effective density of electrons u which 
slides along the channel opening a new view on possi- 
bilities of superconductivity in such materials. 

The question about sliding in such a channel is rather 
nontrivial being linked with fundamental results of dy- 
namical systems [6,7] which we briefly discuss. In fact 
in a local approximation of small charge oscillations the 
forces depend linearly on displacements corresponding 
to a string of particles, linked locally by linear springs 
and placed in a periodic potential. The density of par- 
ticles or charges corresponds effectively to a fixed ro- 
tation number in a dynamical symplectic map which 
describes the recurrent positions of particles in a static 
configuration with a minimum of energy. For linear 
springs this map is exactly reduced to the Chirikov 
standard map [6]. This model in known also as the 
Frenkel-Kontorova model which detailed description is 
given in [8]. At small channel deformation a or small 
amplitude of periodic potential K the particles can 



Preprint submitted to Physica B 



11 October 2011 



Fig. 1. 

freely slide in the periodic potential that corresponds 
to the regime of invariant Kolmogorov-Arnold-Moser 
(KAM) curves which rotation number determines the 
particle density. In this regime the spectrum of small 
excitations is characterized by a phonon spectrum with 
the dispersion relation uj = Csk where k is dimension- 
less wave vector and Cs is dimensionless sound velocity. 
Above a certain critical strength of deformation the 
KAM curve at a given v is destroyed being replaced 
by an invariant Cantor set known as cantori [9] . In this 
regime the excitations above the ground state have a 
gap oj^ = (cek)^ + and the chain becomes pinned 
by the potential. The gap A is proportional to the Lya- 
punov exponent of dynamical orbits on such cantori set 
[7,9] . The transition between sliding and pinned phases 
is known as the Aubry transition [8]. In the pinned 
phase the Aubry theorem guarantees that at fixed f 
there is a unique ground state which static equilibrium 
configuration corresponds to the cantori set with posi- 
tions of particles forming a devil's staircase. However, 
from the physical view point this ground state is rather 
hard to reach since in its vicinity there are exponen- 
tially many equilibrium configurations which energy is 
exponentially close to the energy of the ground state. 
Numerical studies show that already for 100 particles 
the energy difference is as small as 10"'^^ in dimensional 
units [10]. This phase was called the dynamical spin 
glciss (or dynamical instanton glass) since such proper- 
ties appear in spin glasses [11] which have random dis- 
ordered on-site energies and interactions. In contrast 
to that these properties of the Aubry phase appear in 
absence of any disorder being of purely dynamical ori- 
gin of cantori in a purely periodic potential. 

The studies of properties of the quantum Frenkel- 
Kontorova model has been started in [12] and further 
significantly advanced in [13]. It was shown [13] that 
quantum fluctuations lead to melting of the pinned 
phase at sufficiently large values of dimensionless 
Planck constant h. This transition is a zero tempera- 
ture T — quantum phase transition. At small ft ^ 1 
an T = the phonon mode is frozen but quantum 
tunneling gives transitions between quasi-degenerate 
equilibrium classical configurations which can be 
viewed as instantons. At small <g 1 the density of 
instantons is small and their interactions are weak. 
When h increases the instanton density grows and 
above a certain critical ftc ~ 1 the quantum melting of 
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pinned phase takes place at zero temperature leading 
to zero gap, appearance of quantum phonon mode 
and quantum chain sliding [13]. The results obtained 
for the Wigner crystal in a periodic potential [2] in 
classical and quantum regimes confirm this qualitative 
picture. At fixed amplitude of periodic potential K the 
classical Wigner crystal is pinned at small charge den- 
sities f < Vci [2] . Indeed, at ^ we have a problem 
of one electron with zero kinetic energy and obviously, 
an electron is pinned by a periodic potential. 



2. Sliding Wigner snake 

The situation is different in the case of snaked 
nanochannel: noninteracting electrons, corresponding 
to the limit v <^ 1, move freely inside the wiggled 
channel and pinning of the Wigner crystal appears 
only above a certain critical charge density v > 1/2- 
An example of sliding and pinned regimes is shown 
in Fig.2. The data clearly show that the sliding phase 
at a = 0.6 has a smooth hull function and sound 
dispersion law for small oscillations of the crystal. In 
contrast, in the pinned regime at a = 1.2 the hull func- 
tion has a form of fractal devil's staircase and gapped 
spectrum of small oscillations. 

More detailed results on dependence of gap A on 
charge density v and deformation a are described in [3] . 
In [3] it is also shown that for moderate deformations 
a < 1 the charge positions in a static configuration are 
described by a symplectic dynamical map 



V = v + 2a (1 — cos v) sin 20 
(j) = (j) + V + sin v cos 2(j> , 



(1) 



where v — Xi — Xi^i, (j) = Xi axe conjugated action- 
phase variables, bar marks their values after iteration. 
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Fig. 3. 

The map is implicit but symplectic (see e.g. [7]). Ex- 
amples of Poincare sections of this map at two values of 
deformation a are shown in Fig. 3. Phase space region 
with scattered points corresponds to chaotic dynamics 
with pinned phase, while the smooth invariant KAM 
curves correspond to the sliding phase. 

In many aspects the properties of the Wigner crys- 
tal in snaked nanochannels are similar to those of the 
Frenkel-Kontorov model [10] and the Wigner crystal 
in a periodic potential [2]: in the pinned phase there 
are exponentially many static configurations being ex- 
ponentially close in energy and corresponding to the 
dynamical glass phase. However, there are also some 
specific features: for rational values of densities u = 
i^m = where m is an integer, the Wigner snake 

can slide freely since a displacement does not modify 
the Coulomb energy of electron interactions. 

In analogy with the results presented in [2], we ex- 
pect that the quantum Wigner crystal shows a zero 
temperature quantum phase transition going from a 
pinned phase at fi < ftc ~ 1 to a sliding phase at ft > he- 
However, a direct demonstration of this fact requires 
further numerical simulations using quantum Monte 
Carlo methods described in [2,12,13]. 



3. Discussion 

In the above Section we considered the Wigner crys- 
tal in a snaked nanochannel without any internal po- 
tential. It is natural to assume that a more realistic 
case of molecular organic conductors, as shown in the 
Little suggestion in Fig. 1, has not only channel defor- 
mation but also a periodic potential inside the channel. 
Thus the case of organic conductors corresponds to a 
case of snaked channel with a periodic potential inside 
it. The combination of results of [2,3] shows that for 
a given deformation and amplitude of the periodic po- 
tential we have the sliding phase in a certain range of 
charge densities v. 
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We suppose that the sliding KAM phase may corre- 
spond to effective superconducting behavior of electron 
transport in organic conductors. Indeed, the pressure 
diagram of organic conductors shown in Fig. 4 from 
[14] shows that superconductivity exists only is a finite 
range of pressure. We assume that pressure gives vari- 
ation of effective charge density inside the molecular 
channels in the Little suggestion in Fig. 1. This leads 
us to the KAM concept of superconductivity of elec- 
trons without attractive forces: the Wigner crystal of 
electrons slides freely inside a snaked molecular crystal 
channel if the charge density is located inside of KAM 
phase defined by (2). Of course, further studies are re- 
quired for development of this concept. In fact, the slid- 
ing KAM phase can be viewed as a superfluid phase of 
electrons. Indeed, we see that in the KAM phase there 
is a spectrum of excitations with a finite sound veloc- 
ity Cs ■ Thus according to the Landau criterion [15] the 
sliding of electrons with velocities v < Cs is superfluid. 
Hence, the transition from the sliding KAM phase to 
the pinned Aubry phase corresponds to the transition 
from superfluid to insulator. In this superfluid liquid 
the charge carries have change "e" and not "2e" as it 
is the case for BCS pairs. May be effect of interactions 
between electrons in parallel snaked channels should 
be taken into account to have 2e-pairs. We note that it 
is known that repulsive interactions can create super- 
fluid phase in disordered Id systems, e.g. in the repul- 
sive Hubbard model with disorder [16] . 

The existence of dynamical spin glass phase with 
pinned Wigner crystal shows that there should be 
very slow relaxation processes corresponding to very 
slow transitions between quasi-degenerate static equi- 
librium conflgurations. In fact the experiments with 
organic conductors show very slow variations of con- 
ductivity which take place of a scale of days. Such 
experimental results have been reported at ECRYS- 
2011 by K.Miyagawa [17] and P.Monceau [18]. Usually 
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it is argued that the glassy phase appears due to im- 
purities. We think that the origin of this phenomenon 
is not related to disorder and impurities, which pres- 
ence should be rather small in organic crystals used 
in experiments [17,18]. In contrast this glassy phase 
appears as a result of dynamical spin glass phase de- 
scribed in [2,12,13,3] which exists in purely periodic 
structures without any impurities and disorder. 

Finally, wc note that in [2] it was proposed to study 
the dynamical spin glass with cold ions in optical lat- 
tices which can model the problem of Wigner crystal in 
a periodic potential. Such experiments with cold ions 
are now under active discussions [19,20] and their ex- 
perimental realization is on the way [21]. 
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Figure captions 

Fig.l (a) A schematic image of the Little sugges- 
tion for electron transport in organic molecules (after 
[4,5]). (b) A schematic image of electron Wigner crys- 
tal with charges Ci (points) sliding in a snaked sinu- 
soidal nanochannel, dashed lines show force directions 
between nearby electrons. 

Fig.2 Hull function s = h{x) (a,b) and phonon spec- 
trum oj(k/N) (c,d) for incommensurate electron den- 
sity V = N/L = 144/233 is shown at a = 0.6 (a,c) and 
a = 1.2 (b,d). Here x gives the positions Si of electrons 
at o = 0. Here a gives the dimensional amplitude of 
sinusoidal channel described by equation y = a sin x. 

Fig. 3 Poincaro section for the dynamical map (1) at 
a = 0.2 (left panel), 0.6 (right panel). 

Fig.4 Schematic phase diagram of the {TMTCF)2X 
family taken from [14] (see details there) . 



